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摘 要 : 水 资源 短缺 已 成 为 人 类 社会 主要 的 资源 环境 问题 ,气候 和 土地 利用 变化 使 得 水 文 要 素 的 演变 呈现 出 复杂 


性 与 不 确定 性 ,探讨 变化 情景 下 水 文 要 素 空间 分 布 特征 对 区 域 经 济 社 会 可 持续 发 展 具 有 重要 的 理论 和 现实 意义 。 
基于 1960 一 2019 年 大 通 河源 区 气象 水 文 数据 , 通 过 模型 模拟 及 情景 分 割 方法 定量 分 析 了 气候 和 土地 利用 变化 情景 
下 水 文 要 素 空 间 分 布 特征 。 结 果 表 明 :(1) SWAT 模 型 经 过 校准 与 验证 后 ,模型 决定 系数 .纳什 系数 和 百分比 偏差 
(PBIAS ) 均 达到 模型 要 求 , 其 在 率 定期 内 分 别 为 0.81% .0.79% 和 -0.8% ,在 验证 期 内 分 别 为 0.81% 0.75% Fil 15.8% , 


表明 模型 在 大 通 河源 区 具有 较 好 的 适用 性 。(2) 大 通 河源 区 各 水 文 要素 空 间 异 质 性 明显 ,不 能 用 单一 的 水 文 要 素 代 
表 整 体 的 空间 分 布 情况 ;降水 .潜在 蒸 散 发 与 土壤 含水 量 随 海拔 的 升 高 而 降低 ,地 表 径 流 与 产 水 量 随 海拔 的 升 高 而 


增加 。(3) 3 种 情景 下 各 水 文 要 素 的 空间 分 布 大 体 一 致 , 产 水 量 的 空间 分 布 受 土地 利用 变化 情景 的 影响 比较 大 。 气 
候 变 化 情景 使 得 实际 燕 散 发 和 土壤 售 水 量 呈 现下 降 趋 势 , 地 表 径 流 和 产 水 量 呈 现 增 加 趋势 ,土地 利用 变化 情景 下 


各 水 文 要 素 的 变化 与 之 相反 。 
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水 是 地 球 各 种 生命 赖 以 生存 与 发 展 的 基础 ,对 
人 类 有 着 广泛 而 深刻 的 影响 “。 随 着 全 球 工业 化 进 
程 的 不 断 加 速 及 人 口 的 激增 ,导致 水 资源 开发 利用 
不 合理 ,水 资源 供需 矛盾 日 益 加 剧 ,因而 水 资源 短 
缺 问 题 已 逐渐 成 为 人 类 经 济 社会 发 展 所 面临 的 主 


部 的 祁连山 地 跨 我 国 甘 肃 青海 两 省 ,是 西北 干旱 地 
区 最 重要 的 水 源 地 及 生态 屏障 ,其 生态 地 位 十 分 
重要 。 

气候 和 土地 利用 变化 对 水 文 要 素 的 影响 呈现 
出 复杂 性 与 不 确定 性 态势 ,许多 国内 外 学 者 围绕 这 


要 资源 环境 问题 ,受到 国际 社会 的 普遍 关注 。 我 国 
是 水 资源 严重 缺乏 的 国家 ,人 均 水 资源 占有 量 还 不 
到 世界 平均 水 平 的 四 分 之 一 ,并 且 时 空 分 布 极 不 
均匀 ”。 进 入 21 世 纪 以 来 ,极端 王 旱 气候 的 频 发 已 
使 水 资源 匮乏 成 为 影响 我 国 经 济 发 展 、 社 会 安定 和 
生态 环境 改善 的 主要 制约 因素 ”“。 据 统计 ,我 国 每 
年 因 缺 水 造成 的 粮食 损失 达 600 多 亿 斤 "。 中 国 西 
北 干 旱 区 土地 面积 约 占 国土 面积 的 24.5%, 由 于 次 
居 亚 欧 大 陆 腹 地 ,平均 年 降水 量 偏 少 ,是 世界 上 自 
然 环 境 最 严酷 的 干旱 区 之 一 。 位 于 青藏 高 原 东 北 
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两 方面 展开 了 一 系列 的 研究 ,一 方面 ,由 于 自然 
和 人 为 因素 导致 的 全 球 气 候 变 暖 已 成 为 一 个 共 
Wo 2013 年 CC 组 织 的 第 五 次 报告 数据 显示 :过 
去 的 100 a, 地 球 表面 温度 大 约 上 升 了 0.85 % ,而 
国内 众多 学 者 的 研究 结果 都 表明 未 来 一 段 时 间 我 
国 气 温 将 会 呈现 出 继续 增 温 的 趋势 中。 气候 变 瞬 
热 必 会 加 快 水 循环 速率 ,对 不 同时 空 尺 度 上 降水 、 
蒸发 .径流 等 气象 水 文 要 素 产 生 极 大 影响 ” ,并 导 
致 极端 气候 水 文 事件 频 发 中 。 男 一 方面 ,剧烈 的 人 
类 活动 (如 水 利 工程 的 修建 ,城市 化 进程 等 ) 改 变 了 
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流域 下 垫 面 条 件 和 水 资源 在 时 空 尺度 上 的 分 配 过 
程 , 对 流域 土壤 湿度 . 蒸 散 发 和 流域 产 汇流 过 程 造 
成 了 显著 影响 。 而 流域 地 表 径 流 ` 地 下 径流 土壤 
含水 量 等 水 文 要 素 直 接受 降水 因素 的 影响 ”“。 因 
此 ,探讨 气候 及 土地 利用 变化 情景 下 流域 各 水 文 要 
素 的 空间 分 布 特征 对 全 面 理解 水 资源 演变 特征 及 
合理 配置 提供 了 新 的 思路 。 目 前 定量 评估 气候 变 
化 和 土地 利用 对 水 文 要 素 影响 的 方法 主要 有 传统 
经 验 统计 分 析 及 水 文 模型 两 方面 。 由 于 SWAT 模 型 
全 面 的 考虑 了 水 文 要 素 变 化 的 自然 和 人 为 因素 , 充 
分 体现 出 了 空间 异 质 性 ,能 较为 完整 的 揭示 流域 水 
文 过 程 的 内 在 机 理 。 田 唱 等 "基于 SWAT 模 型 评估 
未 来 气候 变化 和 土地 利用 变化 对 汉江 流域 径流 的 
共同 影响 ,结果 表明 气候 和 土地 利用 变化 的 又 加 使 
径流 呈 增 加 趋势 ,气候 变化 较 土 地 利用 对 径流 的 影 
响 更 加 显著 。 任 才 等 "通过 对 叶 尔 羌 河 源 区 径流 
变化 的 影响 因素 人 研究 认为 气候 变化 是 影响 径流 变 
化 的 主要 因素 。 霍 军 军 等 “通过 SWAT 模 型 对 拉萨 
河流 域 径流 进行 研究 发 现 气 候 变化 与 土地 利用 对 
径流 影响 占 比分 别 约 为 82.95% 和 17.05%。 

大 通 河 是 黄河 的 二 级 支流 ,地 处 祁连山 水 源 涵 
养 重要 生态 功能 区 ,是 青海 省 东北 部 和 甘肃 省 西部 
重要 的 水 源 地 ,在 青海 .甘肃 两 省 社会 经 济 发 展 中 
具有 十 分 重要 的 地 位 。 在 全 球 变化 背景 下 ,大 通 河 
源 区 陆地 水 循环 过 程 正 发 挥 着 快速 变化 ,水 文 过 程 
呈现 出 不 稳定 性 和 脆弱 性 的 态势 。 近 年 来 有 许多 
学 者 基于 常规 统计 分 析 方 法 对 大 通 河流 域 水 资源 
状况 展开 了 一 系列 的 研究 ,例如 , 黄 维 东 等 根据 
实测 洪水 资料 分 析 了 梯级 水 电 开 发 对 流域 洪水 过 
程 的 影响 ,认为 水 电站 数量 的 增加 使 得 洪水 过 程 变 
得 极为 频繁 。 董 军 等 他 认 为 1994 年 开始 大 通 河 上 
游 径 流量 缓慢 增加 ,下 游 则 快速 下 降 。 王 大 超 等 " 
采用 回归 分 析 法 对 大 通 河 干流 进行 研究 发 现 
1956—2015 年 大 通 河中 上 游 年 径流 量 以 0.0539x 
10 ma 的 速率 增加 ,部 将 站 则 以 0.0945x10 m°- a” 
的 速率 减少 。 但 这 些 研究 多 围绕 单一 的 径流 要 素 
采用 常规 统计 方法 从 时 间 尺 度 开 展 , 对 大 通 河源 区 
各 水 文 要 素 的 空间 分 布 特征 尚 不 清楚 ,无 法 全 面体 
现 水 文 要 素 空间 异 质 性 。 明 晰 气候 及 土地 利用 变 
化 情景 下 流域 各 水 文 要 素 的 空间 分 布 特征 对 全 面 
理解 水 资源 演变 特征 及 合理 配置 提供 了 新 的 思 
路 。 基 于 此 ,本 研究 以 大 通 河源 区 为 研究 区 ,基于 


1960 一 2019 年 气象 水 文 资 料 建立 SWAT 模 型 ,定量 
分 析 气 候 与 土地 利用 变化 情景 下 水 文 要 素 空 间 分 
布 特征 ,研究 结果 为 区 域 水 资源 管理 及 经 济 社会 可 
持续 发 展 提供 科学 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

大 通 河流 域 位 于 青海 省 东北 部 , 属 黄河 二 级 支 
流 ,河流 发 源 于 祁连山 南 苑 天 峻 县 木 里 境内 , 沿 祁 
连 山 南 莲 向 东 流 经 青海 .甘肃 两 省 的 9 个 县 ,最 后 由 
民 和 县 境内 汇 入 涅 水 河 , 全 长 554 km。 本 研究 所 指 
的 大 通 河 源 区 为 青石 嘴 水 文 站 以 上 部 分 ,流域 内 最 
低 海拔 为 2884 m, 最 高 海拔 为 5067 m ,海拔 较 高 , 落 
差 大 、 地 形 地 貌 复 杂 多 样 (图 1)。 气 候 主 要 受 西 风 
带 、 高 原 季 风 及 夏季 风 的 影响 , 据 统计 资料 显示 
1960—2019 年 多 年 平均 降水 量 为 413.34 mm ,多 年 
平均 气温 为 -0.59 %。 土 地 利用 类 型 以 草地 为 主 , 占 
总 面积 的 50% 以 上 。 

大 通 河 年 均 径 流 约 为 28.95x10: mi, 水 资源 相对 
丰富 ,流域 内 目前 已 建成 跨 流域 调 水 工程 3 处 ,分 别 
为 1994 年 的 “ 引 大 入 秦 ”, 据 统计 年 引水 量 达 4.43x 
10° my ,建成 于 2003 年 的 “引流 济 金 ”工程 的 年 引水 
量 也 达到 了 0.4x10:m;2016 年 “ 引 大 济 涅 ”工程 顺利 
完工 并 通 水 ,年 引水 量 3.5x10: mi, 此 外 地 方 政府 还 
在 考虑 将 流域 内 的 水 资源 调 往 黑河 ,金昌 的 西 大 河 
及 青海 湖 等 地 区 ”7。 

1.2 数据 来 源 

本 研究 所 用 到 的 数据 主要 分 为 栅 格 数据 和 统 
计数 据 两 大 类 ,机 格 数据 有 DEM .土地 利用 土壤 类 
型 等 ,这 些 也 都 是 模型 运行 的 关键 数据 ;统计 数据 
包括 气象 数据 与 水 文 数据 等 ,气象 与 水 文 数据 即 可 
以 用 于 模型 运行 与 验证 ,同时 又 是 历史 气象 水 文 规 
律 分 析 的 基础 数据 。 数 据 详 细 信 息 如 表 1 所 示 。 
1.3 研究 方法 
1.3.1 情景 设置 ”水 文 要 素 的 变化 同时 受到 气候 与 
土地 利用 变化 的 影响 ,为 了 定量 分 析 气 候 与 土地 利 
用 变化 对 水 文 要 素 的 影响 ,根据 1962 一 2019 年 气象 
数据 ,为 保证 气象 数据 时 间 序 列 与 土地 利用 数据 相 
对 应 ,以 1980 年 为 界 将 1962 一 1980 年 、1981 一 2019 
年 气象 数据 与 1980 年 .2020 年 土地 利用 数据 相互 交 
又 ,设计 出 S0.S1.S2 和 S3 共 4 种 情景 ( 表 2) 来 量化 气 
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图 1 大 通 河源 区 区 位 及 子 流域 划分 
Fig. 1 Datong River source region location and subbasin division 
#1 SWAT 模型 所 需 数据 及 来 源 
Tab.1 Data and sources for the SWAT model 
数据 名 称 数据 来 源 数据 描述 
DEM 地 理 空间 数据 云 (http://www.gscloud.cn) 数字 高 程 模型 ,30 mx30 m 
土地 利用 资源 环境 科学 与 数据 中 心 (https://www.resdc.cn) 1980 年 .2020 年 土地 利用 类 型 图 ,1000 mx1000 m 
土壤 数据 联合 国 粮食 及 农业 组 织 (FAO) 土地 利用 空间 分 布 数据 ,1:1 000 000 
气象 数据 中 国 气象 数据 网 (http://data.cma.cn/) 1960—2019 年 祁 连 门 源 .野牛 沟 . 刚 察 , 托 勒 气象 站 气温 E 
水 H RE .风速 .湿度 观测 数据 ,逐日 
水 文 数据 青海 省 水 文 与 水 资源 局 1960 一 2016 年 青石 嘴 水 文 站 径流 观测 数据 ,逐日 


表 2 气候 与 土地 利用 变化 情景 设置 


Tab.2 The setting of scenario for climate and land use changes 


情景 情景 类 别 气候 数据 时 段 土地 利用 数据 
S0 基础 情景 1962 一 1980 年 1980 年 
S1 气候 变化 情景 1981 一 2019 年 1980 年 
S2 土地 利用 变化 情景 1962 一 1980 年 2020 年 
S3 协同 变化 情景 1981—2019 年 2020 年 


候 与 土地 利用 对 水 文 要 素 空间 分 布 特征 的 影响 ”。 
1.3.2 SWAT 模 型 SWAT (Soil and Water Assess- 
ment Tool) 模 型 具有 解决 在 大 尺度 复杂 的 土壤 类 型 、 
土地 利用 对 水 文 水 资源 影响 的 问题 能 力 ””。 模 型 
TE ArcGIS 10.2 平 台 上 基于 流域 的 地 形 、 土 地 利用 类 
型 及 土壤 等 基础 数据 可 以 划分 成 若干 水 文 响 应 单 
元 (HRUs) ,然后 在 每 个 水 文 响应 单元 内 部 根据 水 量 
平衡 原理 计算 相应 的 值 ,最 后 求 得 流域 总 出 水 口 断 
面 的 流量 。 模 型 建立 后 ,需要 经 过 适用 性 评价 才 可 
以 用 来 模拟 大 通 河源 区 径流 情况 所。 本 研究 在 梳 


理 前 人 评价 体系 的 基础 上 ?23 ,选择 确定 系数 (局 )、 
Nash-Suttcliffe( 纳 什 效 率 系数 , 记 作 NS) 和 百分比 偏 
差 (PBIAS)3 个 指标 作为 本 模型 的 评价 标准 。 


2 结果 与 分 析 


2.1 SWAT 模 型 的 率 定 与 验证 

设置 1960 一 1961 年 为 预 热 期 ,1962 一 1999 年 为 
率 定 期 ,2000 一 2016 年 为 验证 期 ,以 1980 年 土地 利 
用 类 型 为 背景 ,结合 研究 区 气温 与 降水 数据 完成 
SWAT 模 型 的 构建 。 在 DEM 集 水 面积 阔 值 为 10000 
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hm 的 情景 下 共生 成 了 39 个 子 流域 。 同 时 ,选取 26 
个 参数 在 swat-cup 中 通过 SUFI-2 算 法 对 模型 进行 
率 定 ,SWAT 模 型 评价 结果 如 表 3 所 示 。 

由 表 3 可 知 ,校准 期 的 结果 为 R°=0.81, NS= 
0.79 , PBIAS=-0.8% ,模拟 值 较 实 测 值 大 ,分 别 为 
50.96 m°- s~! FI 50.56 m’ -+ s~'; 验证 期 的 结果 为 R= 
0.81 , NS=0.75 , PBIAS=15.8% ,模拟 值 小 于 实测 值 ， 
分 别 为 49.97 m-s F 59.33 mes, WRENS 
数 带 回 到 SWAT 模 型 中 重新 运行 便 可 得 到 模拟 的 径 
流量 ,模拟 与 观测 的 径流 量 对 比如 图 2 所 示 ,观测 值 
与 模拟 值 具 有 一 致 的 趋势 ,整体 模拟 效果 较为 理 
想 。 综 上 可 知 ,模型 经 过 率 定 和 验证 后 在 大 通 河源 
区 具有 很 好 的 适用 性 。 

2.2 大 通 河源 区 水 文 要 素 空间 分 布 特征 

根据 率 定 后 的 SWAT 模 型 得 到 39 个 子 流 域 水 
文 要 素 空 间 分 布 特征 (图 3) ,年 降水 量 空间 分 布 在 
467.52~614.41 mm, 低 值 区 主要 分 布 在 流域 西北 部 
和 中 部 部 分 地 区 ,呈现 出 从 东南 向 西北 减少 的 趋 
势 。 这 主要 由 于 流域 东部 地 区 较 中 西部 地 区 更 容 
We 年 实际 蒸 散 

空间 分 布 特征 整体 上 与 降水 空间 分 布 类 似 , 说 


ar s 间 分 布 存在 一 定 的 关联 关系 ,但 在 
流域 中 部 地 区 空间 异 质 性 较为 明显 ,主要 受 其 潜在 
蒸 散发 空间 分 布 的 影响 。 土 壤 含水 量 的 空间 分 布 


没有 明显 的 规律 ,其 高 值 区 主要 分 布 于 流域 东部 和 
中 部 部 分 地 区 ,其 主要 原因 是 流域 东部 地 区 降水 


多 ,中 部 部 分 地 区 实际 蒸 散发 较 少 ,而 地 表 径 流量 


的 空间 分 布 与 产 流量 的 空间 分 布 高 度 趋同 。 综 上 
所 述 ,流域 各 水 文 要 素 空 间 异 质 性 明显 ,不 能 用 单 


一 的 水 文 要 素 代 表 整 体 分 布 情况 。 

从 不 同 海拔 高 程 带 划 分 的 结果 来 看 ( 表 4) BE 
水 .实际 蒸 散发 .土壤 含水 量 均 随 海拔 升 高 而 降低 ， 
而 地 表 径 流 与 产 流 量 随 海拔 升 高 而 增加 ,潜在 蒸 散 
发 的 分 布 无 明显 规律 。 流 域 各 海拔 高 程 带 的 干旱 
指数 相差 不 大 ,空间 分 布 不 明显 ;但 其 蒸发 系数 随 
海拔 升 高 而 降低 ,说 明 在 降水 量 充 足 的 情况 下 28 
散发 能 力 限 制 了 实际 蒸 散发 的 发 生 。 
2.3 气候 与 土地 利用 变化 情景 下 水 文 要 素 空 间 变 
化 分 析 
23.1 气候 变化 情景 下 径流 空间 变化 分 析 气候 变 
化 期 (1981 一 2019 年 ) 与 气候 基准 期 (1962 一 1980 
年 ) 相 比 , 气 温 、 降 水 分 别 增加 了 0.92 % 和 38.90 
mm。 根 据 情 景 So 和 情景 S1 模拟 结果 可 知 ,气候 变 
化 导致 实际 燕 散 发 和 土壤 含水 量 分 别 下 降 8.37 mm 
和 3.53 mm, 使 得 地 表 径 流量 和 产 水 量 分 别 增加 
1.53 mm 和 4.44 mm。 由 气候 变化 情景 下 各 水 文 要 
素 模 拟 结果 与 降水 量 的 空间 分 布 特征 可 知 (图 4)， 
实际 燕 散 发 与 基准 期 相 比 ,除了 19 号 21 号 子 流域 
外 ,其 他 子 流域 均 略 呈现 下 降 趋势 ;土壤 含水 量 从 
原来 的 9.42~251.54 mm 下降 到 8.72~246.94 mm, 仅 
仅 9 号 和 11 号 子 流 域 略 有 增加 ;地 表 径 流 大 的 子 流 
域 主 要 位 于 流域 西北 部 ,尤其 是 1 号 2 号 .3 号 、4 


表 3 大 通 河源 区 率 定期 与 验证 期 模拟 结果 


Tab.3 Simulation results of Datong River headwater rate in periodic and validation period 


时 间 段 观测 平均 值 /ms 站 模拟 平均 值 /(mi…s') NS R PBIAS/% 
率 定期 50.56 50.96 0.79 0.81 -0.8 
验证 期 59.33 49.97 0.75 0.81 15.8 
_ 350 一 观测 值 --- 模 拟 值 
T 300 i 
“=, 250 i 
E 200 | i TE 1 
a ahl 
i | | 
€ AU AMMA Mi AVA AWWA 
a of 
Ee Ne 


日 期 /年 -月 


图 2 大 通 河源 区 模拟 径流 与 实测 径流 对 比 


Fig. 2 Comparison of simulated runoff and measured runoff in the source region of Datong River 
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号 .8 号 .10 号 .14 号 及 16 号 子 流域 ,其 多 年 平均 地 
表 径 流 超 过 203.23 mm, 这 主要 与 其 实际 葵 散 发 低 
和 坡度 大 有 关系 ; 产 水 量 空间 分 布 特征 与 地 表 径 流 
整体 上 具有 高 度 的 一 致 性 。 

总 体 上 ,大 通 河源 区 多 年 平均 地 表 径 流量 最 大 
的 子 流域 与 产 水 量 大 .实际 蒸 散发 小 和 土壤 含水 量 
小 的 子 流域 均 分 布 在 流域 西北 部 ,多 年 平均 地 表 径 
流量 和 产 水 量 大 的 子 流域 分 布 受 流域 西北 部 实际 
蒸 散发 小 .坡度 大 和 地 势 陡 峭 的 影响 ;而 土壤 含水 
量 的 空间 分 布 特征 与 地 表 径 流 和 产 水 量 的 分 布 并 
不 一 致 ,说 明 不 能 用 单一 的 水 文 要 素 代表 整体 分 布 
情况 。 
2.3.2 土地 利用 变化 情景 下 径流 空间 变化 分 析 大 
通 河源 区 主要 土地 利用 类 型 为 草地 , 占 总 面积 的 
50% 以 上 ;其 次 是 裸 地 , 占 总 面积 的 16% 以 上 。 耕 


3 期 赵 美 亮 等 :气候 及 土地 利用 变化 对 大 通 河 源 区 水 文 要 素 空 间 分 布 的 影响 385 
99°E 100°E 101°E 99°E 100°E 101°E 99°E 100°E 101°E 
= Z 多 
名 Z% Z 多 Z 
& 8 <a & 
降水 量 /mm 实际 蒸 散发 /nm 土壤 含水 量 /mm 
四 467.52~487.93 =206.21~231.87 8.95~31.08 
m487.94~515.93 加 231.88~323.37 =31.09~51.17 
515.94~536.72 m323.38~402.29 =51.18~73.61 
= 536.73~602.03 wm 402.30~456.84 加 73.62~91.43 
m602.04~614.41 m456.85~869.66 0 60 km mn 91.44~248.34 
= 
99°E 100°E 101°E 99°E 100°E 101°E 99°E 100°E 101°E 
99°E 100°E 101°E 99°E 100°E 101°E 
Z 
30 = 
Ki 2 
地 表 径 流 /mm 
=18.10~20.37 器 67.58~113.16 
=20.38~77.47 wa 13.17~172.73 
77.48~131.14 四 172.74~220.23 
wm 131.15~203.13 四 220.24~245.52 
四 203.14~303.10 四 245.53~303.46 
99°E 100°E 101°E 99°E 100°E 101°E 
图 3 子 流域 尺度 上 水 文 要 素 空间 分 布 
Fig.3 Spatial distribution of hydrological elements at sub-basin scale 
表 4 大 通 河源 区 不 同 海拔 带 水 文 要 素 分 布 特征 
Tab.4 Distribution characteristics of hydrologic elements at different elevations in the source region of Datong River 
海拔 /m 降水 /mm ” 洪 在 蒸 散 发 /mm 实际 蒸 散 发 /nm 土壤 含水 量 /mm 地表 径 流 /mm 产 流量 /mm ”干旱 指数 蒸发 系数 
2884~3427 568.01 737.38 378.89 63.29 84.86 180.41 1.30 0.67 
3427~3704 546.87 769.00 368.72 61.73 116.44 188.36 1.41 0.67 
3704~3935 535.94 745.05 336.06 52.99 138.11 204.11 1.39 0.63 
3935~4177 536.65 744.37 334.63 52.53 136.26 203.01 1.39 0.62 
4177~5067 535.81 749.81 332.99 51.92 134.83 203.12 1.40 0.62 


地 冰川 \ 居 民 地 和 其 他 建设 用 地 只 占 很 小 的 比例 ， 
不 足 1%。 与 1980 年 相 比 ,2020 年 耕地 、 居 民 区 ,水 
域 和 其 他 建设 用 地 有 所 增加 ( 表 5) ,其 中 ,其 他 建设 
用 地 增加 最 明显 ,增加 了 42.67 km’, 主要 由 草地 转 
化 而 来 ;永久 性 冰川 雪 地 减少 最 明显 ,减少 了 49.08 
km’, 主要 向 裸 地 进行 转移 。 

土地 利用 变化 对 水 文 要 素 的 影响 趋势 与 气候 
变化 情景 相反 ,使 得 实际 蒸 散发 与 土壤 含水 量 分 别 
增加 3.42 mm 和 0.08 mm, 地 表 径 流量 与 产 水 量 分 别 
减少 4.41 mm 和 2.42 mm( 图 5)。 根 据 基 础 情景 对 比 
可 知 ,地 表 径 流量 与 产 水 量 发 生变 化 的 区 域 主要 有 
4 号 .6 号 7 号 及 36 号 子 流域 ;根据 子 流域 主要 土地 
利用 类 型 变化 情况 可 知 ( 表 6),4 号 .6 号 .7 号 子 流 
域 地 表 径 流 及 产 水 量 均 呈 明显 的 下 降 趋 势 , 除 受 薰 
散发 增加 的 影响 外 , 主要 与 裸 地 转换 成 草地 , 而 草 
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99°E 100°E 101°E 99°E 100°E 101°E 
名 名 
ka JZ a JZ 
A Gd on 
实际 蒸 散 发 /mm 土壤 含水 量 /mm 
国 202.01~229.88 号 8.72~30.21 
国 229.89~321.94 国 30.22~S0.28 
321.95~393.35 æ 50.29~72.51 
m 393.36~453.38 æ 72.52~90.10 
m 453.39~869.57 0 60 km mm 90.11~246.79 60 km 
区 Z 
2 Z A Z 
3 > 
地 表 径 流 /mm 
局 8.23~22.03 国 67.73~113.13 
国 22.04~77.46 = 113.14~176.39 
ma77.47~131.39 ma 176.40~224.82 
= 131.40~203.22 j "h m 224.83~250.43 
= 203.23~303.82 0 60 km m 250.44~304.17 
下 = 
99°E 100°E 101°E 99°E 100°E 101°E 
图 4 气候 变化 情景 下 子 流域 尺度 上 水 文 要 素 空 间 分 布 
Fig.4 Spatial distribution of hydrological elements at sub-watershed scale under climate change scenario 
表 5 1980 一 2020 年 土地 利用 类 型 变化 转移 矩阵 
Tab.5 Transfer Matrix of land use type change from 1980 to 2020 /km 
2020 年 
ae a ; 其 他 建 ma f 3 Seas 二 
草地 耕地 居民 区 林地 裸 地 设 用 地 湿地 水 域 ”永久 性 冰川 雪 地 ”总计 
草地 3701.71 4.5 1.62 66.3 52.4 38.92 49.52 64.57 0.07 3979.61 
耕地 0.7 8.87 0.07 0.45 = - - 0.07 - 10.16 
居民 区 0.13 0.06 0.95 0.04 - - - 0.03 = 1.21 
林地 65.32 1.33 0.9 991.19 8.83 0.17 0.64 24.82 = 1093.2 
裸 地 116.7 = = 4.33 1228.33 3.29 9.57 1.67 0.37 1364.26 
其 他 建 213 - - 011 = 1.84 - - - 4.08 
19804F UE FA 
湿地 60.55 0.01 0.54 15.84 0.92 1027.49 7.11 = 1112.46 
水 域 22.8 0.07 = 5.63 1.48 1.61 6.7 254.88 0.03 293.2 
永久 性 = = = = 8.3 = - - 4.81 13.11 
冰川 雪 
地 
总 计 3970.04 14.83 3.55 1068.59 1315.18 46.75 1093.92 353.15 5.28 7871.29 


SM Sls Beat es DDT AK oD HA EHER; AE 36 SF 
流域 菊 散 发 增加 量 较 大 ,但 其 地 表 径流 和 产 水 量 呈 
增加 趋势 ,由 此 可 见 , 草 地 相对 于 林地 更 容易 产生 
地 表 径 流 。 此 外 ,19 号 子 流域 虽然 蒸 散发 量 下 降 
28.53 mm, 但 其 地 表 径 流 和 产 水 量变 化 不 大 ,这 由 于 


该 流域 主要 土地 利用 类 型 为 水 域 有 关 , 占 该 子 流域 
面积 的 77.96%。 

2.3.3 协同 变化 情景 下 径流 空间 变化 分 析 ”根据 协 
同 变化 情景 与 基础 情景 对 比 可 知 ( 图 6) ,实际 蒸 散 
发 .土壤 含水 量 和 地 表 径 流量 分 别 下 降 了 4.87 mm, 
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表 6 主要 子 流域 土地 利用 类 型 变化 情况 及 水 文 要 素 变 化 
Tab. 6 Changes of land use types and hydrological factors in main sub-basins 
子 流 域 编码 ”实际 蒸 散发 /nm ”土壤 含水 量 /mm ”地 表 径流 /mm ” 产 水 量 /mm ”草地 转 林 地 /km ”林地 转 草地 /km 裸 地 转 草 地 /km? 
4 26.93 11.89 -41.43 -27.81 0.00 0.00 14.61 
6 24.04 11.30 -45.95 -28.59 0.00 0.00 2.53 
7 38.39 19.11 -76.58 -48.23 0.00 0.00 637 
36 68.87 5.31 7.78 9.62 0.17 0.24 0.00 
99°E 100°E 101°E 99°E 100°E 101°E 
名 Z 
% lz | z 
en Gd 
实际 燕 散 发 /nm 土壤 含水 量 /mm 
加 214.84~243.17 12.98~32.78 
国 243.18~342.02 国 32.79~S2.97 
国 342.03~414.40 æ 52.98~76.10 
mm 414.41~532.26 m 76.11~95.58 
m 532.27~841.31 0 60 km 95.59~203.23 60 km 
-= 上 三 王政 
99°E 100°E 101°E 99°E 100°E 101°E 
99°F 100°E 101°E 99°E 100°E 101°E 
A A 
Oo oo 
cm on 


地 表 径 流 /mm 
加 0.40~18.99 
国 19.00~89.47 

= 89.48~142.11 
æ 142.12~201.65 
m 201.66~301.63 


99°E 


100°E 101°E 


38°N 


38°N 


产 水 量 /mm 
国 66.79~124.20 

æ 124.21~171.15 
m 171.16~221.23 
m 221.24~253.11 
m 253.12~301.99 


99°E 


100°E 101°E 


图 5 土地 利用 变化 情景 下 子 流域 尺度 上 水 文 要 素 空间 分 布 


Fig. 5 Spatial distribution of hydrological elements at sub-watershed scale under land use change scenario 


3.46 mm Fil 2.92 mm, 产 水 量 增加 1.86 mm。 实 际 蒸 
散发 和 土壤 含水 量 最 大 的 区 域 均 位 于 19 号 子 流域 ， 
尤其 是 实际 蒸 散发 远 超 降水 量 , 达 到 了 840.32 mm, 
土壤 含水 量 也 达到 了 198.48 mm, 主要 是 因为 该 子 
流域 土地 利用 类 型 为 水 域 , 占 比 77.96%。 该 情景 

地 表 径 流量 、 产 水 量 与 气候 变化 情景 不 同 的 是 4 号 
子 流域 呈现 减少 的 趋势 ,分别 为 196.97 mm 和 
239.56 mm, 主要 受 土地 利用 变化 的 影响 ,该 子 流域 
在 1980 一 2020 年 裸 地 转 为 草地 和 水 域 转 为 草地 的 
面积 比较 大 ,分 别 为 14.61 km 和 1.64 km ,而 草地 面 
积 的 增加 使 得 该 区 域 涵养 水 源 的 能 力 增 强 。 此 外 
协同 变化 情景 下 其 他 水 文 要 素 的 空间 分 布 变化 不 
大 ,主要 是 因为 研究 区 土地 利用 变化 不 大 导致 其 对 


水 文 要 素 的 影响 比较 小 。 
3 讨论 


水 文 要 素 变 化 的 影响 因素 主要 体现 在 气候 变 
化 和 土地 利用 变化 两 方面 。 气 候 变 化 对 水 文 要 素 
的 影响 主要 是 由 气温 和 降水 共同 引起 ,降水 是 河川 
径流 补给 的 重要 组 成 部 分 ,而 气温 通常 认为 是 通过 
引起 冰雪 消融 对 高 寒 地 区 的 河流 起 补给 作用 ;气温 
对 径流 的 调控 作用 还 表现 在 通过 影响 蒸 散 发 作用 
进而 减少 径流 量 ”。 由 于 研究 区 位 于 河流 源 区 ,地 
势 较 高 ,气温 在 气候 变化 情景 下 仪 上 升 至 -0.27 C, 
对 实际 蒸 散发 的 变化 影响 较 小 ,因此 气温 引起 径流 
减少 的 因素 可 以 忽略 ,但 对 于 其 季节 性 的 差异 仍 有 
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99°E 100°E 101°E 99°E 100°E 101°E 
Z| Z 
% Z & |Z 
o0 OO 
isa] Cn 
实际 蒸 散 发 /mm 
加 202.00~236.48 œ 12.24~30.04 
= 236.49~336.18 = 30.05~52.89 
= 336.19~405.22 æ 52.90~72.30 
m 405.23~523.95 æ 72.31~92.03 
m 523.96~840.32 0 60 km mm 92.04~198.48 0 60 km 
—— | 
99°E 100°E 101°E 
99°E 100°E 101°E 
z A 
2 Zz & Z 
A A 


地 表 径 流 /mm 
-0.30~22.08 
=22.09~68.29 
™68.30~138.98 
= 138.99~203.07 
ma 203.08~303.83 


99°E 100°E 101°E 


国 67.32~128.05 
国 128.06~176.39 
国 176.40~219.33 


m219.34~249.99 
m250.00~304.17 


99°E 100°E 


图 6 协同 变化 情景 下 子 流域 尺度 上 水 文 要 素 空 间 分 布 


Fig.6 Spatial distribution of hydrological factors at Sub-catchment scale under co-change scenario 


竺 后 续 研 究 。 土 地 利用 是 仅 次 于 气候 变化 影响 水 
文 要 素 的 另 一 因素 ,其 主要 表现 在 不 同 土地 覆 被 对 
水 文 过 程 各 要 素 的 影响 。 而 土地 利用 变化 不 仅 受 
人 类 活动 的 影响 ,还 受气 候 变 化 所 带 来 覆 被 改变 的 
控制 ,其 影响 机 理 和 强度 也 不 尽 相同 。 

在 未 来 研究 中 需要 结合 不 同时 间 尺 度 下 气候 
变化 详细 分 析 不 同类 型 土地 利用 变化 对 流域 水 文 
要 素 的 影响 ,以 便 更 好 地 为 大 通 河源 区 水 资源 合理 
利用 提供 支撑 。 


4 结论 

本 研究 基于 1960—2019 年 大 通 河源 区 气象 水 
文 数据 ,通过 模型 模拟 及 情景 分 割 方法 定量 分 析 了 
气候 和 土地 利用 变化 情景 下 水 文 要 素 空 间 分 布 特 
征 。 得 到 以 下 结论 : 

(1) SWAT 模 型 在 大 通 河源 区 的 月 径流 模拟 效 
果 较 好 ,其 校准 期 的 结果 为 R=0.81, NS=0.79, 
PBIAS=-0.8% ;了 验证 期 的 结果 为 RR=0.81,NS=0.75， 
PBIAS=15.8% ,模型 适用 于 流域 气候 和 土地 利用 变 
化 对 各 水 文 要 素 的 影响 研究 。 


(2) 大 通 河源 区 多 年 平均 地 表 径 流量 与 产 水 量 
随 海拔 的 升 高 而 增加 ,其 最 大 的 子 流域 主要 分 布 在 
流域 的 西北 部 ,流域 各 水 文 要 素 空 间 异 质 性 明显 ， 
不 能 用 单一 的 水 文 要 素 代表 整体 分 布 情况 。 

(3) 大 通 河源 区 土地 利用 类 型 以 草地 为 主 , 气 
候 变 化 和 人 类 活动 共同 导致 1980 一 2020 年 土地 利 
用 类 型 变化 ,主要 表现 在 建设 用 地 和 水 域 的 增加 及 
裸 地 和 冰川 的 减少 ,但 由 于 研究 区 草地 和 林地 禾 盖 
面积 大 ,其 水 源 涵养 能 力 较 强 。 大 通 河源 区 水 文 要 
素 的 变化 主要 受气 候 因 素 的 影响 ,气候 与 土地 利用 
变化 情景 对 各 水 文 要 素 的 影响 呈现 相反 的 趋势 。 
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Effects of climate and land use change on the spatial distribution of 
hydrological factors in the source region of Datong River 
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Abstract: Water shortage has become a major resource and environmental challenge worldwide. Climate and 
land use change have made the evolution of current hydrological factors complex and uncertain. Exploring the 
spatial distribution characteristics of hydrological factors under dynamic scenarios is of vital theoretical and 
practical significance for the sustainable development of regional economy and society. The meteorological and 
hydrological data of the Datong River source region from 1960 to 2019 were used in this study to quantitatively 
analyze the spatial distribution characteristics of hydrological elements under climate and land use change 
scenarios based on model simulation and scenario segmentation. The results showed that: (1) After calibration 
and verification of the SWAT model, the coefficient of determination, Nash coefficient, and percentage bias 
(PBIAS) all met the model requirements of 0.81%, 0.79%, and —0.8% in the rate period, and 0.81%, 0.75%, and 
15.8% in the validation period, respectively, which indicated that the model had good applicability in the 
headwaters of the Chase River. (2) Obvious spatial heterogeneity of hydrological elements was detected in the 
headwater area of Datong River, and a single hydrological element could not represent the overall spatial 
distribution. Precipitation, potential evapotranspiration, and soil water content decreased with the increase in 
altitude, while surface runoff and water yield increased with the increase in altitude. (3) The spatial distribution of 
hydrological factors under the three scenarios were generally consistent, while the spatial distribution of water 
yield was greatly affected by the land use change. Under the climate change scenario, the actual evapotranspiration 
and soil water content showed a downward trend, while the surface runoff and water yield showed an upward 
trend. Under the land use change scenario, the changes of hydrological elements were contrary to these 
observations. 


Keywords: hydrological elements; spatial distribution; SWAT model; Datong River 


